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We investigate a model for the dynamics of a solid object, which moves over a randomly 
vibrating sohd surface and is subject to a constant external force. The dry friction between 
the two soHds is modeled phenomenologically as being proportional to the sign of the object's 
velocity relative to the surface, and therefore shows a discontinuity at zero velocity. Using 
a path integral approach, we derive analytical expressions for the transition probability of 
the object's velocity and the stationary distribution of the work done on the object due to 
the external force. From the latter distribution, we also derive a fluctuation relation for the 
mechanical work fluctuations, which incorporates the effect of the dry friction. 



I. INTRODUCTION 

This paper is a sequel to a previous paper p[], which studied the motion of a solid object (e.g., 
a coin) moving with friction over a horizontal, uniformly rough plate, which is vibrated laterally 
by external Gaussian white noise. Here, we extend this study by adding a constant external field, 
so that we effectively consider the object's motion on an inclined plate with the added influence of 
gravity. The evolution of the velocity v{t) of the object relative to the plate is taken to be given 
by the Langevin equation 

mv{t) + av{t) + a{v{t))AF = F + ^{t), (1) 

where F is the external force associated with gravity and ^(t) is an external noise term associated 
with the random acceleration of the plate. As in p!] , we assume that this noise is a Gaussian white 
noise characterized by 

m) = o, {mm) = ^mt-t'i (2) 

where the brackets indicate an average over the noise, and D is the noise strength. 

The friction between the moving object and the plate is phenomenologically described in Eq. ([T]) 
by two friction coefficients: a (fluid-like) dynamical friction coefficient —av, with strength a, which 
is the solid analogue of the viscous (Stokes) friction in fluids, and a static (dry) solid-solid friction 
—a{v)AF, where Aj? is the strength of the dry friction and the sign function cr(v) is defined as 
equal to +1, 0, — 1 for f > 0, t; = 0, u < 0, respectively. The latter force is often referred to as 
Coulomb friction,^ and was studied in the context of Brownian motion by de Gennes |3J. The main 
feature of this force, which is singular at f = 0, is that it can lead to the sticking of the object on 
the plate when v = and the total external force F + ^ in Eq. ([T]) is smaller in magnitude than 
the "contact" force Ap. When v ^ or + > Ap, then the dry friction simply adds to the 
viscous friction force, and so contributes to the sliding motion of the object (slip motion). 

Since the object is under the influence of two types of external forces, namely the random 
vibrations of the underlying plate, ^(i), and the external force, F, it can be considered to be in a 
nonequilibrium state. In the following, we will investigate the statistical properties of this state by 



^ The history of solid-solid or dry friction goes back a long time; see, e.g., P] and references therein. 
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obtaining analytic expressions for the transition probability or propagator f{v,t\vo,tQ), which gives 
the probability that the object has a velocity v at time t starting from an initial velocity vq at time 
to < t. The propagator is obtained using path integral techniques in the low-noise limit D ^ 0, and 
yields, in the limit t — to ^ oo, a, stationary velocity distribution Ps{v) [1], which characterizes the 
nonequilibrium steady state (NESS) of the object, reached when it moves for a long time compared 
with its intrinsic relaxation time. In this limit, we also obtain the large deviation approximation 
of the probability distribution P{Wr) of the mechanical work Wr performed by the external force 
F on the object during the time interval [t, t + r]. This work Wr is defined as 



and is obviously a random variable, since it is a functional of the random velocity v{s). From 
the large deviation approximation of P(WV), we then obtain a so-called fluctuation relation for 
the work fluctuations, which characterizes the combined effect of the dry friction and the fluid- like 
friction in the NESS. 

The motion of two solids over each other with friction is a ubiquitous problem in nature, of 
great practical and theoretical importance, which occurs, e.g., in geology, physics, engineering and 
biology. Compared with viscous or fluid friction, which is relatively well understood from the 
microscopic point of view, dry friction is quite a complicated process, which may be modeled at 
different scales using different physical scenarios (e.g., plastic deformations, breaking of microscopic 
bonds, etc. [2]). The term —AFa{v) should obviously be considered only as a phenomenological 
or effective model of dry friction. It is a very simplified model when compared to naturally- or 
technologically-occurring situations, but one which incorporates some of the basic features of dry 
friction, in particular, the possibility of alternating stick-slip motion. From the point of view of 
nonequilibrium systems, the model of Eq. ([T]) can also be thought of as a simple phenomenological 
model with nonlinear friction, for which the properties of the work fluctuations in a NESS can be 
investigated theoretically and experimentally. We comment on possible experimental realizations 
of this model in the concluding section in the paper. 



The mathematical technique that we will use to obtain the transition probability /(u, t|fo, to) of 
Eq. ([T]) was already described in [Ij. Here we briefly repeat the basis of this technique, and recall 
a few properties of the deterministic equation describing the motion of the object in the absence 
of noise, i.e., .^(t) = 0. These properties will be useful for understanding all the steps involved in 
the derivation of analytical expressions for f{v,t\vo,to). 



In order to obtain the transition probability f{v,t\vQ,to), we use a path integral approach in 
which f{v,t\vQ,to) is formally expressed as an integral over all paths v{s) leading from an initial 




't+T 



(3) 



II. PRELIMINARIES 



A. Path integral formalism 




where 




(5) 
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is the effective potential asssociated with the friction forces and the constant external force present 
in the Langevin equation, Tm = m/a is the characteristic inertial time scale, A = Ap/m is 
the threshold acceleration associated with the dry friction force, and a = F/m is the acceleration 
associated with the external force. In terms of U{v), the stationary distribution of Eq. Q, obtained 
in the asymptotic limit t — )• oo, takes a Boltzmann-like form 

Ps{v) = Ne-^^^^\ (6) 

where 7 = rm? / D and is a normalization constant [5j. As for f(v,t\vQ,tQ), it is expressed as [6] 

where A is a functional of v{s), 



A[v{s),v{s)]= / C{ij{s),v{s))ds, (8) 

J to 

which will be referred to as the action associated with the path v{s), and C is an Onsager-Machlup 
Lagrangian [71 18] given by 

ms),v{s))=^^[v{s) + U'{v{s))]\ (9) 

In Eq. ([7]), the integral J T>v{s) denotes an integral over all paths v{s) from (^0,^0) to {v,t). The 
Jacobian J[w(s)] originates from the transformation ^{t) — t- v{t) in Eq. Q and is a functional of 
v{s) PIO]: 

J[^;(5)]=e^4^"(^(^))'^^ (10) 

Following our previous work [T], we will evaluate the path integral ([7]) in the limit 7 — >• 00, 
using the saddlepoint approximation. In this limit, the dominant contribution to the path integral 
is due to a particular path v*{s) which maximizes the exponent in Eq. ([T]), or equivalently, which 
minimizes the action A. As a result, v*{s) must satisfy the Euler-Lagrange (E-L) equation 

0. (11) 



ds dv* dv* 

The path v*{s) corresponds to the path with the highest probability among all paths connecting 
(fo, to) and (v, t), and is called, for this reason, the optimal path or simply the most probable path 
connecting (fo,to) and {v,t)?' Once this path is found, the transition probability is approximated 
as 

f{v,t\voM^e-'^^^''''''\ (12) 

where we neglect the contribution of the Jacobian as 7 — )• 00.^ 

The next sections present the solution of the E-L equation in full detail. When choosing initial 
conditions for this equation, it is important to note that, without a, the dynamics described by 



^ In the following, we indicate optimal paths v*{s) with the superscript * and drop the argument s whenever possible 

without confusion, so that then v*{s) is just abbreviated to v* . 
^ The Jacobian may be significant for finite 7 and has then to be taken into account, e.g., when comparing the 

saddlepoint approximation Eq. ( 12 1 with numerical results. 
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Eq. Q exhibits a statistical forward/backward symmetry, in the sense that the probabihty of 
observing the object's velocity v at time t for a given initial velocity vq at to is identical to that 
of observing —v at t for the initial velocity —vq at to ilj- The presence of the constant external 
acceleration a breaks this symmetry, and so one has to distinguish between positive and negative 
Vq for a given a. However, the symmetry is recovered upon changing both vq — )• —vq and a — )• —a, 
so that in the following, we will only consider a > without loss of generality, and will distinguish 
between the two cases f o > and vq < 0. 



B. Deterministic motion 

Before we treat the random motion given by the Langevin equation Q, it is instructive to 
consider the deterministic version of that equation obtained by setting ^(t) = 0: 

v{t) = — -v{t)-a{v{t))A + a. (13) 

Two types of motion arise when solving this equation: 

(a) Stick: When a < A, for both v{t) > and v(t) < 0, the combined acceleration on the right 



hand side of Eq. (13) always acts opposite to the motion, so that the v{t) = axis is an attractor 
of the dynamics (cf. Fig. [ijl^a)). This implies that the deterministic path effectively follows the 
V = axis, even though v{t) = is formally not a solution of Eq. (13). Physically, this behavior 
is expected: if the vibrating plate is sufficiently rough and the external force sufficiently weak, the 
object should get stuck, i.e., reach a zero velocity for any initial velocity vq. 

(b) Slip: When a > A on the other hand, the combined acceleration due to the dry friction 



and the external force, —a{v)A + a in Eq. (13), is always positive. This means that in this case 
the object eventually moves away from the v = axis and approaches then a constant non-zero 
velocity attractor Vss (cf. Fig. [T]|^b)), which is given by 

Vss = {a-A)Tm. (14) 

Importantly, the approach to the zero velocity attractor of the motion is in a finite time, as in the 
force free a = case (cf. Fig. [l]^a) and [1]), while the non-zero attractor is approached with an 
exponential decay over an infinite time. 

In the following, we will encounter many situations similar to case (a) above where the motion 
of the object is well defined for t; > or v < 0, but not at v = because of the singularity of ct{v). 
The procedure that we shall adopt to treat such situations is to invoke continuity: physically, the 
velocity of the object cannot jump in time, so that v{s) must be a continuous function of s. In 
practice, this means that the motion of the object should be determined by considering its equation 
of motion in the regions f > and f < separately, and by matching at the origin the solutions 
obtained, if needed. Figure [l]^b) shows an example of a solution of Eq. ( 13 ) obtained from this 



procedure when vq < and a > A. In this case, the solution must physically cross the v = 



axis in a continuous way to reach the positive steady-state velocity Vss, given by Eq. (14). As the 
solution goes from the lower to the upper half of the {v, s)-plane, it changes slope discontinuously 
because of the sign change of the dry friction term. Such a discontinuity in the slope of v{s) will 
often be encountered in what follows. 



III. OPTIMAL PATHS IN THE VELOCITY-TIME PLANE 



We present in this section the solution of the E-L equation (11) associated with the Lagrangian 
Following the previous section, we write the E-L equation for the optimal path v*{s) as two 
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FIG. 1: Two examples of the object's velocity v{t) described by Eq. (13). (a) The case a < A. For both 
positive (solid curve) and negative (solid dashed curve) initial velocities, the object gets stuck after a finite 
time tst, i.e., v{t) = for all t>tst- At t^t, i) is discontinuous due to the a{v) singularity in Eq. ( 13 ). (b) The 
case a > A. For both positive (solid curve) and negative (solid dashed curve) initial velocities, the object 
approaches the constant velocity Vss, Eq. (14), (thin dashed line) for times larger than the relaxation time 
Tjn- When crossing the v = axis, the slope v{t) is discontinuous. Parameter values: = 1.0, A = 0.7. 
(a) a = 0.3. (b) a = 1.2. 



equations: 

V* - ^v* - —{A-a) = v*>0 (15a) 
V* - ^v* + —{A + a) = v*<0, (15b) 

which are supplemented by the requirement that v{t) be continuous across the boundary u = 0. 
To treat this condition explicitly, we distinguish next between optimal paths which always remain 
on one one half plane and optimal paths which cross the singular w* = axis to span the two half 
planes: 

(a) If the boundary conditions of the E-L equation are such that the optimal path remains 
entirely either on the upper [v* > 0) or on the lower half {v* < 0) of the velocity-time plane, 
respectively, one finds the optimal paths simply by solving either Eq. ( |15a ) or Eq. (15b). The two 
solutions are 

4(s) = S±e^/"™ + C7±e-^/"™ + (a ^ A)r^, (16) 

where the subscripts + and — refer to the upper and lower half plane, respectively. The prefactors 
B± and C± are determined from the boundary conditions 

vl{to)=vo, vl{t) = v (17) 

and are straightforwardly given by 

^ e*/"'"(z; + (A ^ a)Tm) - e*"/"-(z;o + (A ^ a)r^) 

^ _ e*/"-(^o + (A ^ a)r^) - e*p/"-(^ + (A ^ a)r^) 

p{t-io)/r™ _ p-{t-to)/T^ ■ ^^^^ 
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Since the basic solutions are always away from the v* = axis, i.e., never have a zero velocity, 
they represent physically a pure slip motion of the object. We denote such optimal paths as direct 
paths on one half plane (cf. Fig. [2^). 

(b) If the boundary conditions are such that the optimal path has to traverse the singularity 
at the V* = axis (e.g., if the initial velocity is positive and the final velocity negative), v* will 



consist of the basic solutions vl_ and , Eq. (16), on the upper and lower half of the (u *, s)-planes. 



respectively. In addition, at the singular v* = axis, where the optimal path crosses from the 
upper to the lower half plane or vice versa, we impose the condition that the optimal path is 
continuous, since on physical grounds the object's velocity should be continuous. 
Two types of behavior at the singularity v = are then possible: 



1. The optimal path crosses the v = axis only at a single point in time (cf. Fig. [2p). Opti- 
mal paths with a single crossing point are denoted as direct crossing paths and correspond 
physically to a pure slip motion of the object just as direct paths. 

2. The optimal path spends a finite time at v* = 0, i.e., it reaches the v* = axis from the 
initial point (t'o,to)i ^-nd stays on this axis for a finite time, before reaching the final point 
{v,t) (cf. Figs.[2|: and[2]i). This behavior is only possible for a < A, because in that case the 



V* = axis is an attractor of the deterministic motion (see the discussion of Eq. ( 13 ) above). 
Since the deterministic dynamics always has zero action, a solution of the deterministic 
equation also has to represent an optimal path, because the optimal paths are those paths 
with a minimal action. Therefore, we have to admit optimal paths that follow the v = 
axis for a finite time. These optimal paths are denoted as indirect paths and correspond 
physically to a stick-slip motion. 

The behavior of these different types of optimal paths is described quantitatively in the next 
two subsections for the cases a < A and a > A, respectively. For this discussion, we first require 
the actions associated with the basic solutions Substituting Eq. (16) into Eqs. ^ and ([s]) 



yields for the action for paths in the upper half plane and lower half, respectively: 

^ [e*/""'((A T a)rm + \v\)- e*»/--((A T a)r^ + |t^o|)]' 
^^""i'^^iJ 2r^ (e2</-™ - e2*o/r^) 

= K±{v,t;vo,to). (20) 

Due to the presence of the constant external acceleration a, the action of the optimal paths fj. are 
different on the upper and lower half planes. 

A. The case a < A 

In this case, direct paths as well as indirect paths appear for > and < (cf. Fig. [2]). 
The direct paths on one half plane will be labeled by a subscript d and can be represented by 
(cf. Fig.[2|i) 

v*As)Cto = ^Us)\Z,, (21) 



where the u^- are given by Eq. (16) and the boundary conditions vo,to and v,t are explicitly 
indicated. The subscripts -|- and — refer here to direct paths (d) on the upper and lower half plane, 
respectively. These direct paths have the associated action 



(22) 
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FIG. 2: Examples of optimal paths v*{s) in the case a < A for vq > 0. The behavior for vq < follows by 
reflection of the paths at the v* = axis, (a) Direct path on the upper (+) half plane, (b) Direct crossing 
path on two half planes, where the crossing time is denoted by . At u* = the curvature changes from a 
positive to a negative value and the derivative v* is discontinuous, (c), (d) Indirect paths: for sufficiently 
large times the optimal path follows the attractor v* = for a finite time. At the crossovers to/from v* — 
at the times ta+ and tb±, respectively, the slope v* is discontinuous. The ± on the times ta and t}, refer to 
the upper and lower half planes, on which the paths are, respectively. Parameter values: Tm = 1-0, A = 0.7, 
a = 0.3. 



which follows immediately from Eqs. (21) and (20). 

The direct crossing paths will be labeled by a subscript dc and can be represented by (cf. Fig. [2}: 

\0,t> 

(23) 



V*dcis) 




to<S<t'' 
t"" < s <t. 



In Eq. (23), the time when the direct crossing path crosses the v* = axis is indicated by t^, 
which must satisfy to < < t and can be found from the principle of minimal action. The action 
of the direct crossing paths is 



(24) 



which follows from Eq. (jSj) with Eqs. (23) and (20). Thus, the crossing time is given by the two 
equations 

d 



dV 



-[A±(0,t^t;o,to) + AT(^'>*;0,t' 



0, 



(25) 
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which each have a unique real root to < < t. Note that throughout this paper the functions 
A±{v, t;vo,to) are adapted to the case at hand by replacing v, t, vq, to hy the appropriate velocities 
and times. At the crossover time the derivative of the optimal path is discontinuous (cf. Figs. [2^ 
and [2)3). 



The E-L equations (15) imply that for a < A the curvature of the optimal path changes from 
a positive to a negative value upon crossing the v* = axis. This implies that the optimal path 
is then always convex or concave on the upper or lower half plane, respectively. We will see later 
that this change in curvature leads to forbidden regions as in the force free case fT], in the sense 
that there exist regions in the {v* , s)-plane that cannot be reached by any direct path from a given 
initial velocity. Final velocities inside the forbidden region can only be reached by indirect paths. 
This is discussed in detail in Sec. IIVI 

Indirect paths consist of three parts: a first part from the initial point (uq, to) to the v* = axis 
at (0, ta±), a second part along the zero axis from (0, ta±) to (0, th±), and a third part from (0, tb±) 
to the final point {v,t) (cf. Figs. [2j: and[2jl). The first part is given either by v"^ or for vq > 
and Vq < 0, respectively, under the boundary conditions (fc^o) and (0,ta±)- The third part from 
V* = axis is given by v"^ or for vq > and vq < 0, respectively, under the boundary conditions 
{0,th±) and {v,t). 

Indirect paths on one half plane (indicated by the subscript id) are thus parametrized by 
(cf. Fig.[2t 



is) 





vUs) 



\0,ta± 
\vo,t() 

\0,tb± 



to < S < ta± 
ta± < S < tb± 
tb± <S < t, 



(26) 



respectively, while indirect crossing paths (indicated by the subscript idc) are parametrized by 
(cf. FigM 



{s) 







\v(,,to 
\V,t 



to< S < ta± 
ta± < S < tb± 
tb± ^ S < t. 



(27) 



Clearly, the times ta±, tf,± have to satisfy the conditions to < ta± < ti,± < t. The actions 
associated with the indirect paths on one half plane and the indirect crossing paths follow from 



Eq. (26) and Eq. (27) with Eq. (20) 



A[vid^v*a\ = A±{0,ta±;vo,to) + A±{v,t;0,tb±) 
Mv*dc^v*ac] = A±{0,ta±;vo,to) + A^{v,t;0,tb±), 



(28) 
(29) 



where the times ta± and tb± are then determined by the principle of minimal action. This is a 
minimization problem subject to the inequality constraint 



to < ta± < tb± < t, 



(30) 



and has two possible solutions depending on the relative positions of the initial and final velocities 
(cf. [U for an analogous discussion). These solutions are as follows: 

(i) Setting the derivatives of A[v*^, v*^] and A[v*^^, v*^J with respect to ta± and tb± equal to zero 
yields 



ta± = to + Tm In 1^1 + 
tb± = t - In ( 1 + 



l^^ol 



(A =F a)T, 
\v\ 



(A ^ a)r„ 



(31) 
(32) 
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Note that the times ta± are just the times at which a deterministic path, described by Eq. (13), 



would relax to a zero velocity starting from {vo,to). Likewise, tb± are the times at which a deter- 
ministic path starting at {v, t) would relax to a zero velocity, moving backward in time. We remark 
that ta± depends only on vq and not on v, while tiy± depends only on v and not on vq. 



The action of the indirect paths can then be calculated explicitly by substituting Eqs. (31 ) and 



(32) into Eq. (28). This yields the simple result 



Wid^v*d] = A±(0,ta±;vo,to) + A±(z;,t;0,tb±) 
= A±(v,t;0,tb±) 
= U{v). 



Likewise, 



A±(0,ta±;vo,io) + AT(^>*;0,tb^) = U{v). 



(33) 



(34) 



The first term, A-i-(0, ta±; ^^o, io), on the right hand side of the first line of Eq. (33) vanishes, because 
it is the action associated with the deterministic relaxation of the object and for such paths the 
Lagrangian and the corresponding actions are identically zero (cf. Eq. ([9])). The third line follows 
upon substituting Eq. (32) into Eq. (20) and noting that A+(t;, t; 0, 4+) only contributes for positive 
V and A_ (v, t; 0, only for negative t;, giving rise to U{v). 

(ii) If the initial and final velocities of the optimal path are such that th± < ta± in confiict with 
Eq. (poj) , the minimum of Eq. ( 28 ) has to be determined under the equality constraint ta± = t6± . 



If the initial and final velocities have opposite signs the minimization then yields Eq. (25) and the 
solution is in this case, i.e., the indirect crossing paths become direct crossing paths. 



B. The case a > A 

In this case, the constant acceleration overcomes permanently the static (dry) friction, so that 
the object will always exhibit a pure slip motion, while stick-slip motion, represented by indirect 



paths, does not exist (cf. Fig. [3j). The attractor of the dynamics is then given by Vss, Eq. (14), 
which is always positive due to our convention a > 0. If i — io is sufficiently larger than the inertial 
relaxation time r^, the optimal path will follow the attractor Vgs before reaching the final velocity. 
One has then to distinguish between three different types of optimal paths depending on the signs 
of initial and final velocities: 

1. For vq > and f > 0, the optimal path is always given by a direct path on the upper half 
plane, parametrized by Eq. (21) (cf. Figs. [3^ andjs]:). 



2. For either vo > and u < or fo < and f > 0, the optimal path is given by a direct 



crossing path, parametrized by Eq. (23), where the crossing time is given as solution of 
Eq. ([25j) (cf. Figs. [3)3 and [3]; 



For vq < and v < 0, the optimal path can either be a direct path on the lower half plane. 



parametrized by Eq. (21) (cf. Fig. [4p,), or it might cross the v* = axis twice (cf. Fig. [4|; 
The latter case is distinct from the previous cases 1 and 2: first, the optimal path crosses the 
axis in order to reach the attractor at positive Vss, and then goes back to the negative final 
velocity. The optimal path is thus parametrized by three parts (the subscript tc denotes 
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FIG. 3: Examples of optimal paths v*{s) in the case a > A and vq > 0. (a) Direct path on one half plane, 
(b) Direct crossing path. Here, the curvature of v* (s) changes sign when crossing the Vss line (thin dashed 
line); cf. Eq. ( [I4| . (c), (d) For sufficiently large times, both the direct path on one half plane and the direct 
crossing path follow the attractor at Vss before reaching their final velocity. Parameter values: Tm = 1-0, 
A = 0.7, a = 1.2. 



here "twice crossing"): 



v*{s)f'^\ to<S<ti 
v+(s)|g|^^ ti < s <t2 



(35) 



Here, the two crossing times ti and t2 are determined by the principle of minimal action. 
The action of the twice crossing path is clearly: 



tc\ 



A_(0, h;vo,to) + A+(0, t2; 0, h) + A_(i;, t; 0, ta), 



(36) 



so that the crossing times are determined by minimization of A[vf^, v^^]. The crossing times 
ti and t2 in this case can only be determined numerically. 
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FIG. 4: Examples of optimal paths in the case a > A if both vq < and v < 0. (a) Direct path on one 
half plane, (b) For sufficiently large times, the optimal path approaches the attractor before reaching the 
final velocity. This leads to a path that crosses the v = axis twice. Parameter values: r,„ = 1.0, A = 0.7, 
a = 1.2. 



IV. ANALYTIC EXPRESSIONS FOR THE TRANSITION PROBABILITY 



In the previous section we have determined the solutions of the E-L equations (15). The 
transition probability follows from these solutions via the saddlepoint approximation Eq. ( |12[ ). 
However, since different types of optimal paths appear as solutions of the E-L equations (e.g., 
direct and indirect paths in the case a < A), we have to determine from these the one minimizing 
the action for given initial and final velocities, (fo,io) and {v,t). As before, we discuss the cases 
a < A and a > A separately. For certain choices of initial and final velocities, it turns out that 
there are forbidden regions in the (u,s)-plane which cannot be reached by direct paths but only 
by indirect paths. 



A. The case a < A 

In this case, direct paths and indirect paths are possible. We can find the realized optimal 
path between given initial and final velocities, either direct or indirect, as the minimum of their 
corresponding actions. 

We consider first direct and indirect paths on one half-plane (cf. Figs. [2^,c and Figs. [3^,c). 
The direct paths have the associated action ^[■(j^,?;^]) while the indirect paths have the associated 
action A[v*^,v*^]. From the condition 

A[v:„v*,]=A[v*„v*,], (37) 



which, using Eqs. (22) and (33), is equivalent to 

U{v) = A±{v,t;vo,to), (38) 

one can derive boundary lines in the (v,,s)-plane for the final velocity v{t), denoted by w'^^t) and 
w~{t),^ such that, for vq > 0, 

A[v*^, < A[ilv2] for 0<v< w+{t), (39) 



The superscripts + and — indicate boundary lines for a positive or a negative initial velocity, respectively. 
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while, for vo < 0, 



A[v*„v*,]<A[i:,v2 



for 



(40) 



Equation (38) leads, with the two equations above, to two quadratic equations for v, which yield 
analytic expressions for w^{t): 



(A =F a)Tm - e 



{t-to)/Trr 



((A =F a)Tm + \vo\) + J{vi + 2(A =F a)Tm\vo\) {e 



2{t~to)/Tm 



1) 



(41) 



This implies that if either < u < w'^it) for > or w~{t) <v<Q for vq < 0, the action of the 
indirect path is lower than that of the direct path, so that the indirect path will be followed. 

Next, we consider direct and indirect crossing paths (cf. Figs. [2}3, d and Figs, [sja, d). In this 
case the condition reads 



(42) 



and leads to two other critical values or boundary lines for v{t), denoted by u^{t) and u (t), such 
that for vq > 0, i.e., for paths crossing from the upper to the lower half plane. 



for 



u'{t) <v <0. 



(43) 



Likewise, for vq < 0, i.e., for paths crossing from the lower to the upper half plane, we have the 
condition 



A[v*idc^V*ac] < Mv*dc,V*dc 



for 



< ^; < u+{t). 



(44) 



Explicit expressions for u {t) can be obtained without actually solving Eq. (42) for v. Indeed, the 



result of the minimization of Eq. (28) shows that, when the initial and final velocities are such 
that tb- = ta+, the indirect crossing path becomes a direct crossing path with th_ = ta+ = . 
The condition tb- = ta+ leads for > to a boundary line u~^{t), while for vq < 0, the condition 
tb+ = ta- leads to a boundary line u~{t). Using Eqs. (31) and (32) both conditions can be solved 
for the final velocity v, leading to the boundary line 



u^{t) = =f(A ±a)T„ 



(A =F a)Tr, 



(A =F a)Tm + \vo\ 



,{t-to)/T„ 



(45) 



It follows from this discussion that for a given initial state {vQ,to), the optimal path is an 
indirect path, either on one half plane or crossing on two half planes, if the final state {v, t) lies in 
the interval u^{t) < v < w~^{t) {vq > 0) or w^{t) < v < u^{t) (vq < 0) with t > ta±, otherwise 
the optimal path is a direct or a direct crossing path (cf. Figs. [5^ and[5}3). The curves u^{t) and 
w^{t) thus represent boundaries in the {v, s)-plane, separating different dynamical behaviors of the 
moving object in terms of direct and direct crossing paths (pure slip motion) and indirect paths 
(stick-slip motion). 

Having determined the boundaries of the different kinds of optimal paths in the {v, s)-plane, 
we now obtain analytical expressions for the transition probability f{v,t\vo,to) by using the sad- 



dlepoint approximation, Eq. (12), i.e., by selecting that optimal path, which minimizes the action. 



For simplicity, we give the results for vq > only. The results for vq < follow similarly. We 
distinguish two cases: 

(a) For t < ta+ no indirect optimal paths can occur, so that, for a final velocity v < 0, the direct 



crossing path v^^, with an action given by Eq. (24), is the optimal path. For v > 0, on the other 
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(a) (b) 




FIG. 5: The boundary lines of the optimal paths for the case a < A indicating a change in the dynamical 
behavior of the object for (a) vq > and (b) vq < 0. If the final velocity v lies either in the interval 
u^{t) < V < w~^{t) {vq > 0) or w~{t) < V < u'[t) (wp < 0), the object follows a stick-shp path since 
the optimal path is either an indirect path on one half plane or an indirect crossing path (the green curve 
indicates an indirect crossing path). Otherwise, the optimal path is either a direct or a direct crossing path 
(blue and red curves, respectively). The boundaries between direct and indirect paths are indicated with 
dashed black curves. Parameter values: = 1.0, A = 0.7, a = 0.2. 



hand, the direct path on the upper half plane, with an action given by Eq. (22), is the optimal 

(46) 



path. By combining these two cases, we thus obtain for the transition probability 
f{v,t\vQ,to) = Ml 



g-7[A+(0,t X ;i)o,to)+A- {v,t;0,t>' )] 



v<0 
v>0, 

in the case t < ta+, where A/i is a normalization factor (cf. Fig. |6^). 

(b) For t > ta+, indirect paths occur and the boundary lines of the optimal paths indicate that 



1. For a final velocity v < u^{t), the direct crossing paths with action given by Eq. (24) are 
optimal; 



2. For u~^{t) < V <0, the indirect crossing with action given by Eq. (34) are optimal; 

3. For < V < w^{t), the indirect paths on the upper half plane are optimal. The actions of 
these paths are given by Eq. (33) and are the same as those of indirect crossing paths; 



4. For V > w~^{t), the direct paths with action given by Eq. (22) are optimal. 
Combining these cases, we then obtain 



f{v,t\vo,to) =M2 < 



(47) 



V > W~^, 



where A/2 is a normalization factor (cf. Fig. |6|3-d). 

Note that the contribution of the indirect paths for < v < is just the stationary distri- 
bution Eq. This implies that we recover the stationary distribution in the asymptotic time 
limit from Eq. (47), 



lim f{v,t\vo,to) =Ps{v), 



(48) 
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FIG. 6: Transition probability f{v,t\vo,0) for a < A and Vq > 0, as given by Eqs. (46) and (471, plotted for 
four different times t. Transitions between the different types of optimal paths are indicated with vertical 
dashed lines. At u = 0, the slope of the transition probability is always discontinuous, (a) For t < ta+, only 
direct crossing and direct paths contribute for v < and v > 0, respectively, in which case f{v,t\vQ,0) is 
given by Eq. (46). (b) For t > ta+^ the distribution is given by Eq. (47) and indirect paths then contribute 
for u'^ < V < ur. At the transition between indirect and direct paths, i.e., at , the slope of the transition 
probability is discontinuous, (c) For larger t, the indirect paths contribute over a larger region of values 
for V. In this figure the contribution of the direct crossing paths is already far in the tail of f{v,t\vo,0), so 
that they do not appear in the figure, (d) For asymptotically large only indirect paths contribute and the 
stationary velocity distribution Ps{v), Eq. is recovered. Parameter values: Tm = 1.0, A = 0.7, a = 0.2, 
7 = 10.0, Wo = 2.0, to — 0, ta+ = 1.6. The values of w+ and w"*" can be read off the figures. 



since then both — >■ oo and u'^ — )• — oo as t — )■ oo, respectively, so that only indirect paths 
contribute to the transition probability in the asymptotic time limit. 

The transition probability f{v,t\vo,to) is plotted in Fig. [6] for four different t values. One of the 
dominant features of the transition probability is the discontinuity in the slope at v = 0, which 
appears for all t values and is a direct consequence of the (t{v) singularity in the equation of motion 
(Q. In addition, the slope of the transition probability is also discontinuous at the boundary line 
w^, where the optimal path switches from an indirect to a direct path, i.e., to a different dynamical 
behavior. 
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B. The case a > A 



It follows from the discussion of the optimal paths given in Sec. |IIIB that, for vq > 0, the 
optimal path is a direct path if the final velocity v is also positive (case 1 in Sec. IIIB), whereas 



the optimal path is a direct crossing path if t> < (case 2 in Sec. IIIB). Therefore, for vq > 0, the 



transition probability is given by Eq. (46) for all t. 



The situation is different for vq < 0. If the final velocity v is positive, then the optimal path 



will just be a direct crossing path (case 2 in Sec. IIIB). However, if < 0, two types of paths 



are possible, namely a direct and a twice crossing path (case 3 in Sec. IIIB). Which one is chosen 



depends on the location of the final velocity and is determined by the principle of minimal action. 
From the condition 



follows then a critical value of v, denoted by Wcr{t), such that 



A[vlv*a]<A[vl,^ 



tcl 



for 



Wcr{t) < V < 0. 



This indicates that 



1. For a final velocity v < Wcr{t), the direct paths with action (22) are optimal; 



2. For Wcr{t) < V < 0, the twice crossing (tc) paths with action (36) are optimal; 



3. For V >0, the direct crossing paths with action (24) are optimal. 



From these results and Eq. (12), we thus obtain the transition probability for vq < 0: 

^ ^-^A-{v,t;vo,to) 

fiv,t\vo,to) = AAs < 



e y < Wcr 

g-7[A-{0,ti;Do,to)+A+(0,t2;0,ti)+A_(t.,t;0,t2)] <v<0 
-7[A_(0,tX;t>o,to)+A+{i>,t;0,t><)] 



> 0, 



where A/"-? is a normalization factor. 



(49) 



(50) 



(51) 



V. THE STATIONARY DISTRIBUTION OF THE MECHANICAL WORK 



In this section, we investigate the properties of the mechanical work, performed on the object 
of mass m due to the external force F = ma in the NESS. This work is given as a functional of 
the object's velocity, Eq. ([3]), and is therefore a random variable. In the following, we calculate the 
probability distribution of the mechanical work Wr performed over a time r > by the external 
force. Since this quantity is expected to scale proportionally to r, we actually consider, from now 
on, the work per unit time or work rate 



w 



Wr/r. 



(52) 



The probability distribution associated with this quantity is denoted by Pr{w), and is calculated 
in the asymptotic time limit r — )• oo. From this calculation, we also derive an important result 
characterizing the NESS, referred to as a fluctuation relation. 
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A. Distribution of the mechanical work 



The path integral representation of Pt{w) is 

P^{w) = [ J[v{s)] e-T^[^W'^W] S {Wr[v{s)] - Tw) Vv{s). (53) 



Similarly as we did for the propagator f{v, t\vo, to), we shall approximate the path integral of Pt{w) 
by its dominant path, which corresponds here to the path that minimizes the action A subject to 
the constraint that the work rate along this path has the value w. The expression of Pt{w) resulting 
from this approximation is expected to be accurate in the regime where both r — )■ oo and 7—7-00, 
i.e., where the integration time r is long and the noise power D of the external vibrations is small. 
To be more precise, we expect, following the theory of large deviations jllj (see also [12] and 
references therein), to obtain an approximation for Pr{w) of the form 

Pr{w) ^ 00,7 ^00, (54) 

where I{w) is a function, called the rate function, which does not depend on 7 or r. Given that 



this approximation is consistent with taking the optimal path of the full path integral of Eq. (53), 
I{w) must be given by 

I{w) = lim -mmA[v{s),v{s)], (55) 

r— )-oo T 

where the minimum of the action, given by Eqs. ([s]) and ([o]), has to be determined among all 
paths v{s) over the time interval [0,t] which give rise to the work value V^V[^^(s)] = wt. No limit 
involving 7 appears in the above expression because 7 is already factored out of the action; see 



Eq. (53). Equivalently, one can express the rate function by 



I^yj) = lim ^^Kl<l^ (56) 

r— >-oo 7" 

if we denote by v^{s) the optimal path that minimizes A subject to the constraint VFr[^^^] = 
As we did for the propagator, the correction that comes from the Jacobian in the path integral is 
neglected as 7 — )• 00. 



To solve the constrained minimization problem of Eq. (55 ), we transform it, following Taniguchi 



and Cohen [13], to an unconstrained minimization problem for the modified action 

Kr[Hs),v{s)] = A[v{s),v{s)] - pWr[v{s)], (57) 

which involves a Lagrange multiplier /3 associated with the constraint Wr[u(s)] = wt. Using 
Eqs. ([s]) and ([3|, the modified action Eq. (57) can be expressed in terms of a modified Lagrangian 



}C{v{s),v{s)) = C{v{s),v{s)) - l3mav{s), (58) 

so that^ 

Kr[v{s),v{s)]= [ }C{ij{s),v{s))ds. (59) 
Jo 



^ We set to = 0, without loss of generality. 
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FIG. 7: Examples of optimal paths that minimize the action K^-, Eq. (57). The different paths are 

distinguished by different initial and final velocities, which can be read off the figure, e.g., the blue curve is 
the optimal path with the initial velocity vq — 2.0 and the final velocity v — 1.0 at the time r = 10, etc. For 
sufficiently large r, the optimal path is close to a constant velocity (dashed horizontal line), which is the 
attractor of the constrained motion described by the modified Lagrangian Eq. (58 1, (a) w = 2.0 and r = 10. 
(b) w = 2.0 and r = 30. The time spent close to the attractor scales with r, while the times that the 
optimal path needs to reach from the initial velocity or to reach the final velocity from -5^, respectively, 
do not. (c) w — 4.0 and r ~ 10. (d) w ~ 4.0 and r = 30. Parameter values: Tm = m = /S. = 1.0, a = 0.5. 



The paths minimizing the modified action Kr are then found by solving the functional equation 
6Kr = or, equivalently, by solving the E-L equation 



d die 
ds dvl 



die 

dvl 







(60) 



for the optimal path v*p{s), where the subscript f3 indicates the dependence of the optimal path 
on the Lagrange multiplier /3. This E-L equation must be solved over a time interval [0, r] with a 
given initial velocity f^(0) = vq and a final velocity v'^i'T') = v. The value of /3 is found by solving 
the constraint equation W^b^] = wt for /3. The Lagrange multiplier /3 is a function of the work 
rate w only, i.e., /3 = (3{w), and the optimal path v^, required to determine the rate function in 



Eq. (56), is given by < = v^^^^. 

Figure [7| shows typical optimal paths that minimize the action K-j- for different initial and 
final conditions as well as for different values of w. The essential property of these paths, which is 
clearly seen in this figure, is that they are attracted to a constant value v^, and deviate from Vw 



18 



only near the initial and final velocities. This behavior of follows because the Euler-Lagrange 



dynamics determined by Eq. (60) has a unique attractor at Vw, whose value depends on w. As a 



result, for sufficiently large r, all optimal paths satisfying the E-L equation (60) converge to the 
attractor starting from any initial velocity vq, then remain near for most of the time, and 
depart from Vw near the end of the time interval [0, r], to join the final velocity v{t) = v. Figure [7] 
also shows that the time that the optimal path needs to reach Vu, from vq as well as to reach v from 
Vw, respectively, do not scale with r. This implies that the time the optimal path spends near the 
attractor scales with r. 

The conclusion that we reach from these observations is that, as r — )• oo, the main contribution 
to the action j4[i;J^, w^] comes increasingly from that part of which follows the attractor ii^^,. For 



the purpose of calculating the rate function I{w), which requires the limit r — >• cxd (cf. Eq. (56)), 
we can therefore assume that the optimal path has the simple form 

vl{s)=v^, (61) 
i.e., is constant for all s, but with the constant depending on the work rate fluctuation w. 



To find the actual value of Vw associated with a given work fiuctuation w, we substitute Eq. (61 ) 
into the definition of the work, Eq. ([s]), to obtain 

Wr[vw] = Tmavw ■ (62) 



Using Eq. (52) therefore leads to the simple result 

w 

Vw = ■ 

ma 



(63) 



The rate function I^w) is then obtained from Eq. (56) by substituting the constant optimal path 



of Eq. (61) with Eq. (63) in the expression of the action, which leads to a constant Lagrangian 
term: 

I{w) =\im-[ C{iw,Vw)ds = C{0,Vw) = ](^^ + a(—)A-a] . (64) 
T^oo T Jo ^ ' ' 4 \maTm \maJ J 

This rate function is plotted in Fig. [sj We can see that I{w) has a discontinuity at li; = 0, which is 
due to the singular nature of the dry friction at u = 0. The jump of I{w) at w = is equal to aA. 

At this point, it is unclear whether the jump of I{w) at w = is a real feature of P{Wt) or only 
an artifact of the large deviation approximation. On physical grounds, we expect this distribution 
to be continuous at w = 0, and it may well be, accordingly, that its large deviation approximation 
must be constructed, as for optimal paths, by considering the regions f > and v < separately, 
and by matching at the origin the results obtained. Settling this issue would involve a study of the 



contribution of the Jacobian correction term, Eq. (10), in the large deviation approximation (54), 
which is beyond the scope of this paper. 

B. Fluctuation Relation for the mechanical work 

Fluctuation relations are mathematical relations for the fiuctuations of thermodynamic quanti- 
ties, such as heat (entropy production) or work in nonequilibrium systems (see, e.g., |14j and ref- 
erences therein). In a NESS, the so called asymptotic or steady state fluctuation theorem (SSFT) 
states that the probability distribution Pt{A) of finding a particular value of a thermodynamic 
quantity A over a time r satisfies a certain symmetry relation of the form (cf. |15H18j ) 



Pr{-A) 



e""^, (65) 
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where c is a constant. The SSFT Eq. (65), also refered to as conventional or Gallavotti-Cohen 
fluctuation theorem |16j . represents a refinement of the second law of thermodynamics in that 
it quantifies the probability of observing fluctuations of thermodynamic quantities around their 
average values in a NESS. 

Using our analytical result for the rate function I{w), Eq. (64), we can derive a fluctuation 



relation for the fluctuations of the mechanical work due to the external force in the stationary 
regime.^ The ratio of positive to negative work fluctuations in this case is just given by the 



difference of the corresponding rate functions (cf. Eq. (54)) 



Prjw) 
Pr{-W) 



^T'y[I{-w)-I(w) 



(66) 



With Eq. (64) one then obtains 



I{-w) - I{w) 



w 



rriTm 



+ aA, 



with the convention w > 0, so that the fluctuation relation is given by 



^r(^^) 
Pr{-W) 



T7| — ^i^+aA 



(67) 



(68) 



This fluctuation relation differs from the SSFT, Eq. (65), by the term TjaA in the exponent, which 



scales linearly with the dry friction coefficient A and leads to a discontinuity at w = of magnitude 
2a A (cf. Fig. Without dry friction, i.e., for A = 0, one recovers the conventional fluctuation 



relation from Eq. ( 68 ) with c given by 



7 



niTr, 



(69) 



This is expected, since for A = our model Eq. ([T]) maps formally onto the paradigmatic nonequi- 
librium particle model, investigated by van Zon and Cohen [T9], for which the conventional work 



^ Importantly, these work fluctuations are macroscopic fluctuations, induced by the mechanical vibrations of the 
plate due to the externally imposed noise. 
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FIG. 9: The difference /(— w) — liw), Eq. (67). The jump at u> = has the magnitude 2aA, i.e., twice the 
jump of I{w) at w = 0. Parameter values: Tm = m = 1.0, A = 1.0, a = 0.5. 



fluctuation theorem was shown to hold. Note, however, that in the latter model the Gaussian noise 
originates from an equilibrium heat bath leading to the constant c = l/fc^T. Note also that, if 
PiyV-j-) were to be continuous because of added corrections coming from the Jacobian term in the 



path integral, then there would be no discontinuity in Eq. (67) or in Fig. |9j as the added term oA 
would simply disappear. 



VI. CONCLUDING REMARKS 



We have investigated a simple phenomenological model Eq. ([T]) for the motion of a solid 
object, which moves over a randomly vibrating solid plate and is subject to a constant 
external force. The most probable or optimal path between an initial state (veto) and a 
final state (f , t) in the (f , s)-plane can be a direct or an indirect path, which correspond, 
respectively, to a slip motion with u 7^ or to a stick-slip motion, where the object is stuck 
on the plate for some time. These two kinds of motion underlie the behavior of many, more 
complicated systems with solid-solid friction. From a mathematical point of view, our model 
also represents one of the simplest models for which singular features of optimal paths can 
be investigated analytically (see also pO]). 



2. To complement our study of the random motion of the solid object, we have investigated 
the fluctuations of the mechanical work performed on the object by the external force. Here, 
the path integral approach allows a straightforward derivation of the rate function for the 
stationary work distribution, which appears to be characterized by a discontinuity at zero 
work due to the singularity of dry friction at = 0. The rate function was used to formulate 
a fluctuation relation for the work fluctuations by taking the ratio of the probability of 
observing positive work values to that of negative work values. The resulting fluctuation 
relation also exhibits a discontinuity at zero work, while behaving linearly with w away from 
zero, similarly to the SSFT. 

3. In addition to the mechanical work due to the constant force, which is given by Eq. (|3]), 
one could also consider the work done on the object by the vibrations of the plate due to 
the externally imposed noise. Over a time r, the amount of work performed due to these 
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vibrations is 

W:;^= r v{s)i{s)ds. (70) 

JO 

An investigation of the fluctuation properties of this quantity requires considerable analytical 
efforts (cf. pT] for the case of Brownian motion), which is beyond the scope of this paper. 

4. In order to be able to compare our analytical expressions with results of numerical simu- 
lations, one has to consider corrections for finite 7, since simulations are performed with 
finite noise strength, while our theory is strictly only valid for 7 00. One correction term 
is given by the contribution of the Jacobian J[f (s)] in the path integral Eq. ([T]), which is 



explicitly given by Eq. (10). This term gives a 7-independent prefactor in the saddlepoint 
approximation Eq. (12) and can only be neglected in the limit 7 — )• 00. Another correction 
term arises from the expansion of the action in the neighborhood of the optimal path v* 
leading to the so-called fluctuation factor F[v*] [lOl [22]. The saddlepoint approximation 
including these two correction terms is then [T] 

f{v,t\voM ^ Jb*]e-^^[^*''^*lF[^*], (71) 



where J[v*] is given by Eq. (10) and the fluctuation factor by 

F[v*] = /■^°'^-7/4d^[^W^+^(^*W)^(^)1p,(,), (72) 

^(0,to) 

with Q.{v) = U" {v)"^ + U' {v)U"' {v) . Both correction terms contain delta-function singularities 
and require careful treatment. In order to avoid ambiguities arising from these singularities, 
one might consider a regularization of the (t{v) term in Eq. ([T]), e.g., by replacing (y{v) by 
tanh(u/e), where e is a small parameter. This way one obtains a nonlinear but not singular 
equation of motion for the object. The singular behavior is then recovered by taking the 
limit e — )• 0.^ 



5. Dry friction depends on the properties of the interface between the two three-dimensional 
solids and is thus an intrinsic two-dimensional effect. It would be interesting to study the 
behavior of two-dimensional generalizations of the one-dimensional model given by Eq. ([T]) , 
possibly with a non-uniform surface roughness, and to compare the results of these models 
with experiments, in order to obtain information about surface properties, such as roughness, 
defects, etc. 

6. Further extensions of the phenomenological model Eq. ([T]) are possible, by adding, e.g., an 
external linear spring force acting on the object. Models of this type, in the absence of 
external vibrations, have been considered as a common type of model in the literature on 
dry friction \12>\ [23]. Another possible extension concerns the statistics of the externally 
imposed vibrations. By considering, e.g., Poissonian shot noise |6] (see also I25j), which 
allows, in principle, a great variety of pulse characteristics absent in Gaussian noise, one 
might obtain new effects in the stick-slip motion of the object. 



This is reminiscent of taking the thermodynamic hmit for a system in thermal equilibrium, which is necessary to 
understand a discontinuous first-order phase transition. In this context, one could also consider the discontinuity 
in the slope of the transition probability /(«, t\vo, to) as analogous to a first-order phase transition in equilibrium. 
The action would then be analogous to the equilibrium free energy. 
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7. The linear dynamic friction —av{t) in Eq. ([T]) is present for a number of solid-solid systems 
for large velocities or for solid-fluid systems (lubricated interfaces). For smaller velocities, a 
different and more complicated nonlinear behavior is often observed, due to the nature of the 
object-solid interface Without going into the meso- or microscopic nature of the velocity 
dependence of these measured dynamical friction coefficients, it would be instructive to 
investigate the phenomenological model Eq. ([T]) , using a simple mathematical representation 
of the measured dynamical friction as a function of the object's velocity. 

In carrying this out, one would study an equation with two different kinds of nonlinearities, 
associated with both dynamical and dry friction. It would be interesting to study the physical 
difference of the system's behavior due to these two nonlinearities. 

8. A recent experimental investigation of the motion of a solid object moving over a rough plate 
subject to externally-imposed, random vibrations and a constant external force (gravity), as 
described by Eq. ([T]), has been performed by Chaudhury et al [271 EH]. Here, the trajectories 
of a small glass prism, which moves over an inclined glass plate in the presence of external 
white Gaussian noise and gravity, are recorded using a high-speed video camera. This 
investigation focuses mainly on the properties of the drift velocity as well as the diffusive 
behavior of the prism, where the diffusion coefficient is estimated from the variance of the 
probability distribution of the prism's displacement. It would be interesting to compare a 
theoretical prediction for the diffusion coefficient of the prism, following from our theory, 
with the experimental results obtained in \27\ [28] . The diffusion coefficient Dp of the prism 
could be evaluated using the Green-Kubo relation (cf. [3]) 



The velocity correlation function {v{0)v{t)) is then obtained from the transition probability 
f{v,t\vo,to) by 



where p{vq) is the stationary velocity distribution Eq. ([6]). 

In [28], the fluctuations of the mechanical work rate, Eq. ([s]), performed by the constant 
external force on the prism have also been determined. In fact, a fluctuation relation that 
increases linearly with the work rate has been measured, which would be in partial agreement 
with our result Eq. (681). However, the jump at zero work rate in the fluctuation relation, 
as predicted by Eq. (68|, was not observed. A more detailed comparison of our theory with 
these experiments is left for the future. 

9. Another interesting experimental test of our theory could be performed in an electrical 
circuit, based on a well-known mapping of a class of mechanical systems on electrical circuits 
[291 [30] . Ill th^-t case, the dry friction can be represented by a diode (cf. [23j). Electric circuits 
are particularly suitable for investigating the properties and physical significance of optimal 
paths (see, e.g., [201 EH |32] ) . 

10. It might be interesting to investigate our model in the context of the recent transition path 
theory |33[ |3^ . Open questions concern the nature of the singularities in the transition paths 
and their connection with optimal paths in the path integral framework. 




(73) 
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